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Abstract 


Given  three  disjoint  n-sets  and  the  family  of  all  weighted  triplets  that 
contain  exactly  one  element  of  each  set,  the  3-  index  assignment  (or  3- 
dimensional  matching)  problem  asks  for  a  minimum-weight  subcollection  of 
triplets  that  covers  exactly  (i.e.,  partitions)  the  union  of  the  three  sets. 
Unlike  the  common  (2-index)  assignment  problem,  the  3-index  problem  is  NP- 
complete.  In  this  paper  we  examine  the  facial  structure  of  the  3-index 
assignment  polytope  (the  convex  hull  of  feasible  solutions  to  the  problem) 
with  the  aid  of  the  intersection  graph  of  the  coefficient  matrix  of  the 
problem's  constraint  set.  In  particular,  we  describe  the  cliques  of  the 
intersection  graph  as  belonging  to  three  distinct  classes,  and  show  that 
cliques  in  two  of  the  three  classes  induce  inequalities  that  define  facets  of 
our  polytope.  Furthermore,  we  give  an  0(n4)  procedure  (note  that  the  number 
of  variables  is  n^)  for  finding  a  facet-defining  cl ique- inequal ity  violated  by 
a  given  noninteger  solution  to  the  linear  programming  relaxation  of  the  3- 
index  assignment  problem,  or  showing  that  no  such  inequality  exists.  v 
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1.  Introduction 

The  (axial)  three-index  assignment  problem ,  to  be  denoted  AP3,  also  known  as 
the  (axial)  three-dimensional  matching  problem,  can  be  stated  as  follows:  given 
three  disjoint  n-sets,  I,  J,  and  K,  and  a  weight  associated  with  each 

ordered  triplet  (i,j,k)  £  I  x  J  x  K,  find  a  minimum-weight  collection  of  n 
disjoint  triplets  (i,j,k)  e  I  x  J  x  K. 

An  alternative  interpretation  of  AP3  is  as  follows.  A  graph  is  complete 

if  all  of  its  nodes  are  pairwise  adjacent.  A  maximal  complete  subgraph  of  a 

graph  is  a  clique.  A  graph  is  k-partite  if  its  nodes  can  be  partitioned  into  k 

subsets  such  that  no  two  nodes  in  the  same  subset  are  joined  by  an  edge.  It 

is  complete  k-partite ,  if  every  node  is  adjacent  to  all  other  nodes  except 

those  in  its  own  subset.  The  complete  k-partite  graph  with  n^  nodes  in  its 

ith  part  (subset)  is  denoted  K  . 

nl»n2* . .  • 

Consider  now  the  complete  tri-partite  graph  Kn^n  with  node  set 
S  =  I  U  J  U  K,  |I|  =  |J|  =  ) K |  =  n.  Figure  1  shows  Kp  n  n  for  n=2  and 
n=3.  K _  n  n  has  3n  nodes  and  n3  cliques,  all  of  which  are  triangles 

n,nTn 

containing  exactly  one  node  from  each  of  the  three  sets  I,J,K.  Let  (i,j,k) 
denote  the  clique  induced  by  the  node  set  {i,j,k}.  If  a  weight  c^  is 
associated  with  each  clique  (i,j,k),  then  AP3  is  the  problem  of  finding  a 
minimum-weight  exact  clique  cover  of  the  nodes  of  Kp^p,  where  an  exact 
clique  cover  is  a  set  of  cliques  that  partitions  the  node  set  S. 


AP3  can  be  stated  as  a  0-1  programming  problem  as  follows: 

max  ^iel^jeJ^keK  CijkXijk 
s,t*  ^ je J^keK  xijk  =  1  V  1eI 
hJktK  Xijk  =  1  v 

y  •  tJ  •  ,  X.  ..  =  1  V  keK 
Mel^JeJ  ljk 

X1  J*k  £  {0,1}  V  i ,  j , k 

where  I,  J  and  K  are  disjoint  sets  with  |I|  =  |J|  =  |K|  =  n.  The  coefficient 
matrix  of  AP3  for  the  case  n=3  is  shown  in  figure  2. 
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Figure  2 


One  interpretation  of  the  constraints  is  in  terms  of  an  nXnXn  cube  in 
three-space  made  up  of  n^  unit  cubes,  each  to  be  assigned  a  value  of  zero  or 
one  so  that  the  following  conditions  hold:  if  the  cube  is  viewed  as  a  set  of 
j,k-planes  stacked  up  in  the  direction  of  the  i-axis,  the  first  set  of 
constraints  requires  that  the  total  value  of  the  variables  in  each  such  plane 
be  exactly  one;  and  a  similar  interpretation  holds  for  i,k-planes  (the  second 
group  of  constraints,  corresponding  to  the  j-axis)  as  well  as  for  i,j-planes 
(the  third  group,  corresponding  to  the  k-axis).  Figure  3  illustrates  this  for 


We  will  denote  by  AP3n  the  (axial)  3-index  assignment  problem  of  order  n 
(i.e.,  defined  for  n-sets),  by  Ap  the  coefficient  matrix  of  its  constraint 
set,  and  by  I  ,  J  ,  Kp  the  3  associated  index  sets.  The  row  and  column  index 
sets  of  An  will  be  denoted  by  Rn  and  Sn  respectively.  Clearly,  |Rn|  =  |Ip|  + 

|Jnl  +  1 Kn t  =  3n  and  sn  =  1 1 n I  x  ! Jn !  x  lKJ  =  "3* 

In  terms  of  Kp  n  n,  Ap  is  the  incidence  matrix  of  nodes  versus  cliques 
(triangles):  it  has  a  row  for  every  node  and  a  column  for  every  clique  of 


n,n,n* 

As  usual,  the  support  of  a  (row  or  column)  vector  is  understood  to  mean 
the  index  set  of  its  nonzero  components.  Each  element  of  S  (that  indexes  a 
column  of  Ap  and  a  clique  of  <n  n  n)  will  also  be  used  to  denote  the  support 
of  the  given  column  of  Ap  and  the  node  set  of  the  given  clique  (triangle)  of 
Kp  n  n.  Thus,  if  as  has  support  (i,j,k)  (i.e.,  if  clique  s  of  Kp  n  n  has  node 
set  (i,j,k})»  we  will  write  s  =  (i,j,k)  or  as  =  aijk,  meaning  that  column  as 
has  ones  in  positions  lei,  jeJ  and  keK. 

AP3  is  a  special  case  of  the  (axial)  3-dimensional  transportation 
problem,  in  which  the  right-hand  sides  of  the  constraints  can  be  any  positive 
integers,  the  sets  I,J,K  are  not  necessarily  equal  in  size,  and  the 
integrality  constraints  are  relaxed.  This  is  in  turn  a  generalization  of  the 
well-known  transportation  problem,  a  special  case  of  which  is  the  simple 
assignment  problem. 

Our  problem  is  called  axial  to  distinguish  it  from  another  3-dimensional 
assignment  problem,  called  planar ,  which  can  be  formulated  as  follows: 

IMX  litl'-jtjikeK  cijkxijk 
s-t.  l)sI  x.jk  ■  1  V  jcJ.keK 

*J.J  *ijk  -  1  * 

IkcK  *ijk  =  1  V 


7 


This  problem  is  a  special  case  of  the  planar  3-dimensional  transportation 
problem.  If  we  view  the  cube  as  a  plane  of  n^  columns  of  unit  cubes,  then  the 
constraints  of  the  planar  problem  require  that  the  sum  of  the  values  of 
variables  in  each  column  be  exactly  one,  in  each  of  the  three  possible 
orientations.  This  is  illustrated  in  Figure  5. 


The  3-dimensional  transportation  problem  (TR3)  in  these  and  other 
formulations  was  first  studied  by  Schell  [20].  The  literature  on  this  problem 
includes  the  references  (2,4,5,9,10,12,13,14,15,18,19,20,211.  The  original 
motivation  for  considering  this  model  was  a  problem  in  the  transportation  of 


> 


3 


goods  of  several  types  from  multiple  sources  to  multiple  destinations. 
Applications  of  AP3  mentioned  in  the  literature  include  the  following 
(Pierskal la  [18,191).  * 

•  In  a  rolling  mill  with  |ll  soaking  pits  (temperature  stabilizing 

baths),  schedule  |K|  ingots  through  the  pits  so  as  to  minimize  idle¬ 
time  for  the  rolling  mill  (the  next  stage  in  the  process).  * 

•  Find  a  minimum  cost  schedule  of  a  set  of  capital  investments  (e.g., 
warehouses  or  plants)  in  different  locations  at  different  times. 

•  Assign  troops  to  locations  over  time  to  maximize  a  measure  of 
capabi 1 ity. 

•  Launch  a  number  of  satellites  in  different  directions  at  different 
altitudes  to  optimize  coverage  or  minimize  cost. 

AP3  is  known  to  be  an  NP-complete  problem  [11].  Obviously,  AP3  is  a 

special  case  of  the  set  partitioning  problem  (SPP): 
max  cx  (1) 

s.t.  Bx  =  e 

X  £  (0.1}q 

where  B  =  (b  )  is  a  matrix  of  zeroes  and  ones  and  e  is  a  vector  of  ones.  The 

set  being  partitioned  in  this  case  is  R  =  I  U  J  U  K,  with  the  rows  of  B 

corresponding  to  the  elements  r  e  R.  The  subsets  s  =  {i,j,k}  from  which  the  ^ 

partition  is  to  be  selected  are  those  subsets  (of  cardinality  three)  which 
contain  one  element  from  each  of  the  sets  I,  J  and  K.  Each  column  bs  of  B  is 
the  incidence  vector  of  one  such  subset  (i.e.,  bs  has  a  one  in  each  row 

correspond i ng  to  an  element  in  the  subset  and  zeroes  elsewhere).  The  set 
packing  relaxation  (SP)  of  (1)  is  the  program: 

max  cx.  (2)  m 


s.t.  Bx  <  e 


For  properties  of  SPP  and  SP  see  the  survey  3]. 

Let  Pj  denote  the  convex  hull  of  feasible  solutions  to  AP3n,  i.e., 

3 

P.  =  convex  £  0,l}n  |A  x  =  e’> 

I  ’  1  n 

The  intersection  graph  GA  =  (V,E)  of  a  0-1  matrix  A  has  a  node  s  for  every 

column  as  of  A,  and  an  edge  (s,t)  for  every  pair  of  columns  as,  such  that 

as«at^0.  The  intersection  graph  Gfl  of  A  is  the  clique-intersect - ^aph  of 

n 

K_  n  n,  i.e.,  Gfl  has  a  node  for  every  clique  (triangle)  of  _  and  an 

n ,  n ,  n  Mp  ..,n,ri 

edge  for  every  pair  of  triangles  that  share  some  node  of  Kn  n.  ^  ,raph  G* 

n  9  ri ,  n  ' 

for  n=2  is  shown  in  figure  6. 
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Although  the  3-4ndex  assignment  problem  has  a  sizeable  literature,  no 
work  has  been  done  on  describing  the  polytope  Pj.  In  this  paper  *e  app'y  the 
tools  of  polyhedral  comb inatorics  to  AP3n  and  obtain  a  partial 
character ization  of  the  facial  structure  of  Pj.  In  particular,  in  section  2 
we  identify  three  classes  of  cliques  of  the  intersection  graph  of  Ap  and  show 
that  they  are  exhaustive.  These  cliques  are  known  to  define  facets  of  the 
polytope 

3 

Pj  -  conv-x  e  <0,1 ,n  |A^x  <  e  ■ , 


» 


> 


> 


the  set  packing  relaxation  of  the  set  partitioning  polytope  Pr.  In  section  3 
we  show  that  two  of  the  3  classes  of  cliques  also  define  facets  of  Pj  (whereas 
the  cliques  in  the  third  class  define  improper  faces),  and  that  these  facets 
are  all  distinct.  Finally,  in  section  4  we  give  an  0(n4)  procedure  for 
detecting  a  clique  inequality  violated  by  some  solution  to  the  linear 
programming  relaxation  of  Pj,  or  showing  that  no  such  inequality  exists. 


2.  The  Cliques  of  GA 

In  this  section  we  identify  all  the  cliques  of  GA,  the  intersection  graph 

of  A. 

For  any  subset  V  r  S  of  the  node  set  of  G^,  we  will  denote  by  <V>  the 
subgraph  induced  by  V.  For  r  £  R,  we  will  denote  by  Sr  the  support  of  row  r 
of  A,  i.e.,  Sr:=  fseS|ars  =  1}. 

Proposition  2.1.  For  each  r-;R,  the  node  set  Sr  'Induces  a  clique  (of 

O 

cardinality  nc)  in  G^. 


11 


Proof.  The  subgraph  <Sr>  is  obviously  complete.  To  see  that  it  is 

maximal,  assume  w.l.o.g.  that  re  I  (an  analogous  reasoning  holds  if  reJ  or 
reK).  Now  let  SeS\Sr  be  arbitrarily  chosen,  and  let  s  =  (iQ,j0,k0).  Since 
s^Sr,  r^iQ;  and  Since  Sr  contains  all  triplets  whose  first  element  is  r, 
there  exists  teSr,  t  =  (r,j,k),  such  that  r^iQ,  j7j'q,  k^kQ.  Hence  Sr  U  {s} 
does  not  induce  a  complete  subgraph  of  GA;  and  since  this  is  true  of  any 
SeS\Sr,  the  subgraph  of  GA  induced  by  Sr  is  maximally  complete,  i.e.,  a 
clique.  Furthermore,  |Sr|  =  n^  for  all  re R . | ( 

The  set  of  cliques  defined  by  Proposition  2.1  will  be  called  class  1  and 
denoted  Q^.  Clearly,  |Qj  =  3n.  In  terms  of  Kp  p  n,  the  clique  of  class  1 
corresponding  to  row  r  of  A  contains  those  nodes  of  the  intersection  graph  GA, 
whose  associated  triangles  in  Kp  n  n  share  node  r  of  Kp  n  p. 

Proposition  2.2.  For  every  seS,  let 
T(s)  =  {t  e  S\{s}|as*at  =  2}. 

Then  the  node  set  {s}  U  T(s)  induces  a  clique  of  size  3n  -  2  in  GA. 

Proof.  Let  s  =  00»j0.k0)»  and  let  tj_,t2  e  T(s)  be  chosen  arbitrarily, 
with  tj^t2.  Since  each  of  t^  and  t2  contains  two  of  the  three  elements 

i0,j0,l<o»  tj  and  t2  must  have  at  least  one  element  in  common.  Hence  the  node 
set  {s >  U  T(s)  induces  a  complete  subgraph  in  GA.  Now  let 
u  s  S\({s>  U  T(s) ) .  Then  the  triplet  u  =  (i,j,k]  contains  at  most  one 
element  of  s.  If  au*as  =  0,  we  are  done.  Assume  now  that  au«as  =  1,  with  i 
=iQ  (a  similar  reasoning  holds  if  j=jQ  or  k=kQ).  Then  j^jQ  and  k^kQ.  Now  by 
definition,  T(s)  contains  some  t  =  (i^*j0,k0)  such  that  i  ^  iQ(=i).  But  then 
au«at  =  0,  i.e.,  (u >  U  (s)  U  T(s)  does  not  define  a  complete  subgraph  of 
Ga.  Since  the  choice  of  u  was  arbitrary,  the  subgraph  defined  by  {s>  U  T(s) 
is  maximal  complete. 


For  each  seS  and  for  each  of  the  three  pairs  of  the  triplet  s  = 
(i0*J0*k0) ’  there  are  n_l  other  triplets  in  S  containing  the  same  pair;  hence 
| T ( s ) |  =  3(n-l),  and  thus  {s}  U  T(s)  has  3n-2  elements. || 

The  set  of  cliques  defined  in  Proposition  2.2  will  be  called  class  2  and 
denoted  Q2.  There  is  exactly  one  clique  of  class  2  for  every  column  of  A,  and 
there  is  no  double  counting;  hence  | Q2 !  =  n^.  In  terms  of  Kn  n  n,  the  clique 
of  class  2  corresponding  to  column  s  =  (iQ,j0,k0)  of  A  contains  the  node  of  GA 
corresponding  to  the  clique  (i0,j0,kQ)  of  Kp^n  n,  along  with  the  3(n-l)  nodes 
of  GA  corresponding  to  those  cliques  of  Kn  n  n  that  share  an  edge  (a  pair  of 
nodes)  with  the  clique  (i0,j0,k0). 

Proposition  2.3.  For  every  ordered  pair  s,t  e  S  such  that  as,at  =  0, 

let  ti,  t2,  t3  e  S  \  {s,tf  be  the  (uniquely  defined)  triplets  such  that 
1  c  tT  t.. 

as-a  1  =  1,  a^a  1  =  2,  i  =  1,2,3. 

Then  the  node  set  {s,t1,t2*t3}  induces  a  (4-)clique  in  GA. 

Proof.  Let  s,t  e  S,  with  as*at  =  0,  and  let  s  =  (is,js,ks),  t  = 

(it,jt,kt).  Then  t1  =  (is,jt,kt),  t2  =  (it,js,kt)  and  t3  =  (it,jt,ks)  are  the 

only  3  triplets  in  S\(s,t}  that  satisfy  the  requirements  of  the  Proposition  , 

s  ti 

i.e.,  they  exist  and  are  unique.  Further,  a  »a  =1  for  i=l,2,3  and 

t.  t. 

a  «a  J  =  2  for  all  i,j  e  {1,2,3};  hence  {s.t^.t^tj}  induces  a  complete 

subgraph  in  GA.  To  see  that  this  subgraph  is  maximal,  note  that  any  triplet 

u  t  S\{s>  that  contains  an  element  of  s,  either  contains  two  elements  of  t 

(and  hence  is  identical  to  one  of  the  triplets  t^,  t2  or  t2),  or  else  contains 

t  • 

at  most  one  element  of  t.  But  then  au»a  1  =  0,  where  ti  £  {t^ , t2,t3)  is  the 
triplet  containing  those  two  elements  of  t  not  contained  in  u  (besides  the 
element  of  s).  Thus  (s,t]_,t2,t3)  induces  a  maximal  complete  subgraph,  hence  a 
4-clique  in  GA. | | 


The  set  of  cliques  described  in  Propositions  2.3  will  be  called  class  3 
and  denoted  Q3.  In  terms  of  Kn>r|fn.  every  class  3  clique  of  GA  is  associated 
with  an  ordered  pair  (s,t)  of  disjoint  triangles  of  Kp  n  n,  and  its  node  set 
contains  (a)  the  node  of  GA  corresponding  to  the  triangle  s,  and  (b)  the  3 
nodes  of  GA  corresponding  to  those  triangles  tj,  t2,  t2  of  Kn>n^n  that  share  1 
node  with  s  and  2  nodes  with  t. 

As  to  the  cardinality  of  Q3,  every  ordered  pair  (s,t)  such  that 
as«at  =  0  gives  rise  to  a  clique  of  class  3.  Since  |S|  =  n^  and  for  every 

O  e  f 

s  e  S  that  are  (n-l)J  indices  t  e  S  such  that  a  *a  =  0,  the  number  of 
ordered  pairs  (s,t)  with  as»at  =  0  is  n^(n-l)^. 

To  determine  the  number  of  cliques  of  class  3  we  also  need  to  know  how 
many  different  ordered  pairs  give  rise  to  the  same  clique.  Let  s  = 

0'St jsiks) »  T  =  ( i t » * *^t ) *  ^1  ~  (^s*^t*^t^*  ^2  ~  (^t*^s*^t^*  ^3  -  ^t’^t’^s^’ 
and  denote  by  C(s,t)  the  node  set  of  the  clique  (of  class  3)  corresponding  to 

the  ordered  pair  (s,t),  i.e.  let  C(s,t):=  {s.tj^^}.  Further,  let 
?!  =  (it,Js,ks),  ?,  =  (ic,  j*,kc) ,  f,  =  (i^j^kj.  Then  we  have 


s’^t’  s 


S,JS*  t; 


Proposition  2.4.  C(s,t)  =  C(t.,E.)  for  i  =  1,2,3. 

Proof.  Consider  the  ordered  pair  (tp?^.  From  the  definitions,  the  4 
triplets  of  the  set  C(t1,?1)  are  (i$,jt,kt)  =  t1#  (is,js,k$)  =  s,  (it,jt,k$) 
=  t3,  and  (it,js,kt)  =  t2;  thus  CftpEj)  =  C(s,t).  By  symmetry,  C(t.,E.) 

=  C(s,t)  for  i  =2,3.  |  | 

Corollary  2.5.  The  number  of  cliques  of  class  3  is  n^(n-l)^/4. 

Proof.  Every  clique  of  class  3  arises  from  4  distinct  ordered  pairs,  and 
the  number  of  the  latter  is  n^(n-l)^.  || 

Next  we  show  that  GA  has  no  other  cliques  than  the  ones  described 
above.  But  first  we  need  a  property  of  Gfl. 


Proposition  2.6.  GA  is  regular  of  degree  3n(n-l)  and  has  ^  n  (n_1) 
edges. 

Proof.  Let  as  be  an  arbitrary  column  of  A.  There  are  (n-1)^  columns  a^ 
of  A  such  that  as»at  =  0,  hence  there  are  n^-l-(n-l)^  =  3n(n-l)  columns  au 

of  A  such  that  as*au  f  0.  Thus  the  degree  of  node  s  in  GA  is  3n(n-l),  and 

by  symmetry  this  is  true  of  all  s  e  S.  Since  the  number  of  edges  of  a  graph 

is  one  half  of  the  sum  of  the  degrees  of  its  nodes,  GA  has 

2  x  n^  x  3n(n-l)  =  jj  n^(n-l)  edges.  || 

Let  Q  denote  the  set  of  cliques  of  GA. 

Theorem  2.7.  The  only  cliques  of  GA  are  those  of  classes  1,  2  and  3; 
i.e.,  Q  =  Qx  U  q2  U  q3. 

Proof.  We  will  use  induction  on  n.  For  n  =  2,  the  statement  is  found  to 
be  true  by  listing  all  the  cliques  of  GA  (see  Fig.  4).  In  fact,  because  of 
symmetry,  it  is  sufficient  to  list  the  cliques  containing  a  given  node,  say  1, 
and  they  are  [1,2, 3, 4},  [1,2, 5, 6},  [1,3, 5, 7}  (of  class  1),  (1,2, 3, 5}, 
(1,2, 4, 6},  (1,3, 4, 7},  (1,5, 6, 7}  (of  class  2),  and  (1,4, 6, 7}  (of  class  3). 

Suppose  now  that  the  statement  is  true  for  n  =  2,...,k  and  let  n  = 

k+1  >  3.  Consider  the  relationship  between  Ga  and  Ga  .  Note  first  that 

Mn-1  Hn 

Kp  n  n  is  obtained  from  Kn_  1 , n- 1 , n- 1  ac^n9  Three  nodes  i*,  j*,  k*,  to 

the  three  sets  I,  J,  K,  respectively,  so  that 

[n  ’  I„-l  U  (1*>.  0n  *  Jn_!  U  (J.)  «•>  Kn  *  Kn_!  U  («.); 

and  (b)  adding  new  edges  (i*,  i)  for  all  i  t  Jn  U  K  ,  (j*,  i)  for  all 

i  c  In  U  Kn,  and  (k*,  i)  for  all  a  e  In  U  Jn. 

This  creates  3n(n-l)+l  new  triangles,  namely: 

(i)  one  new  triangle  (i*,j*,k*)  that  shares  no  node  with  any  triangle  of 

*n-l,n-l,n-l’ 

(ii)  3(n-l)  new  triangles  of  the  form  U,j+,k*)  for  all  i  e  I  ,, 
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(i*,4,k*)  for  all  a  e  J  ,,  and  (i*,j*,a)  for  all  i  e  K  and 

n— i  n- I 

(iii)  3(n-l)^  new  triangles  of  the  form  (p,q,k*)  for  all 

(p,q)  e  In  l  X  Jn_j,  (p, j*,q)  for  all  (p,q)  e  If1_1  x  K^,  and 

(i*,p,q)  for  all  (p,q)  e  Jnl  x  K^. 

In  terms  of  the  coefficient  matrix  of  AP3,  An  is  obtained  from  Af|_1  by 

adding  3  new  rows,  one  to  each  of  the  sets  Ip_lt  Jp_lt  with  zero  entries 

in  the  positions  indexed  by  Sp_^;  and  adding  3n(n-l)+l  new  columns  (of 

dimension  | Rp j  =  |Rp_jJ+3),  with  supports  corresponding  to  the  new  triangles 

of  Kp  n  n  described  under  (iii)  above. 

It  then  follows  that  Ga  is  obtained  from  Ga  by  adding  3n(n-l)+l  new 

Hn-1 

nodes  corresponding  to  the  triangles  described  under  (i),  (ii),  (iii),  and  a 

new  edge  for  every  pair  of  nodes  (s,t)  of  Ga  such  that  (a)  at  most  one  of  s 

rtn 

and  t  is  a  node  of  G»  ,  and  (b)  the  two  triangles  of  K_  _  _  corresponding  to 

Mn-1  nfnfn 

s  and  t  have  at  least  one  node  (of  Kp  p  p)  in  common.  It  also  follows  that 

Ga  is  an  induced  subgraph  of  Gfl  . 

Hn-1  Hn 

Consider  now  the  node  set  C  of  an  arbitrary  clique  of  Gfl  .  The 

Hn 

restriction  C  of  C  to  Ga  induces  a  clique  of  Gfl  ,  hence  by  the  induction 

M  Hn-1 

hypothesis  <C>  belongs  to  one  of  the  three  classes  described  in  Propositions 

2.1,  2.2,  2.3.  With  an  argument  analogous  to  the  one  used  in  the  proof  of  the 

corresponding  Proposition  to  show  that  a  complete  subgraph  of  the  given  class 

is  maximal,  we  will  show  for  each  of  the  3  cases  that  <C>  belongs  to  the  same 

class  as  <C>.  Since  <C>  was  chosen  arbitrarily,  this  will  prove  that  the 

three  classes  of  cliques  are  exhaustive  for  Ga  =  Ga  ,  thus  completing  the 

*n  Hk+1 

induction. 

Suppose  first  that  <C>  e  Q^.  Then  C  is  of  the  form 

C  =  Sp_j.  =  (j  e  S n _ ^ f a j  =  1}  for  some  i  e  Rn  W.l.o.g.  we  may  assume 

that  i  e  I  ,  (an  analogous  reasoning  holds  if  i  e  J  , ,  or  i  e  K  .). 

n-i  J  3  n-I  n-1' 


Now  suppose  <C>  i  Q^.  Then  <C>  has  a  node  s  =  li0,j0,kQ}  such  that  iQ  i 
i.  Since  ^  contains  all  the  triplets  of  Sp_^  whose  first  element  is  i, 
it  certainly  contains  some  t  =  (i,j,k)  such  that  j  f  jQ  and  k  +  kQ.  But  then 
as»at  =  0  (where  as  and  a^  are  columns  of  Ap),  hence  we  obtain  the 
contradiction  that  <C>  is  not  a  clique. 

Suppose  next  that  <C>  e  Q^.  Then  C  is  of  the  form  {s}  u  T(s),  where 

s  e  Sp  x  and  T(s)  =  {t  e  Sn  ^\{s}|as»at  =  2}  (as  and  at  being  columns  of 

Ap_^).  Let  s  =  (iQ,j0,k0).  Now  suppose  <C>  t  Qp.  Then  <C>  has  a  node  u  = 
(i,j,k),  such  that  a5«au  <  1  (where  as  and  au  are  columns  of  A  ).  Since  <C> 
is  a  clique,  as»au  f  0,  hence  as«au  =  1.  W.l.o.g.,  assume  that  the  common 

element  of  the  triplets  s  and  u  is  i'q  (a  similar  reasoning  holds  if  the  common 

element  is  jQ  or  kQ).  Then  j  f  j  and  k  j-  kQ.  By  definition,  T(s)  contains 

some  t  =  (i*,  jQ,  kQ)  such  that  i*  ^  iQ(=i),  ami  t  e  T(s)  implies  t  £  C. 
But  then  au»at  =  0,  contradicting  the  assumption  that  <C>  is  a  clique. 

Finally,  suppose  <C>  e  Q^.  Then  C  is  of  the  form  {s.tj.tg.tj},  where 
s,  tj_,  t2,  t3  e  sn  i  are  distinct  and  such  that 


e  t.  .  t. 

as«a  1  =  1,  ax«a  1  =  2, 


i  =  1,2,3, 


t-  z  ■ 

for  some  t  e  Sp_^  such  that  as»at  =  0.  (Here  as,  at  and  all  a  1  are 

columns  of  Ap_1) .  Let  s  =  (is,js,k$),  t  =  (it,jt,kt),  and  tj_  =  (is,jt,kt),  t2 

=  (it,j$,kt),  t3  =  (it,jt,ks).  Now  suppose  <C>  l  Q^.  Then  <C>  has  a  node 

u  £  S  \S  ,,  say  u  =  (i,j,k),  such  that  au«as  >  0,  au*  a  7>  0,  i  =  1,2,3 
n  n_i  t. 

(here  au,  as  and  all  a  1  are  columns  of  Ap).  Since  u  l  Sp_^,  at  least  one 

of  the  elements  i,  j,  k  is  not  contained  in  any  of  the  4  triplets  s,  tj,  t2, 

t3*  W.l.o.g.  assume  k  is  such  an  element.  Then  each  of  s,  t^,  t2,  t3  must 
contain  either  i  or  j.  There  are  four  cases:  i  =  i$,  j  =  j$;  i  =  it,  j  =  jt; 

i  =  i„.  .i  =  i*:  and  i  =  n.  .i  =  .i_.  In  each  case  at  least  one  of  the  four 


triplets  s,  t^,  t^,  tj  does  not  contain  either  i  or  j,  which  contradicts  the 
assumption  that  <C>  is  a  clique. 

Thus  every  clique  of  G«  =  G*  belongs  to  either  Q, ,  or  Q?,  or  Q-,,  and 

Hn  Hk+1  1  c  J 

the  induction  is  complete.  || 


3.  Facets  of  Induced  by  Cliques  of  GA. 

If  C  is  the  vertex  set  of  a  clique  of  GA<  then  obviously  every  x  e  p 

satisfies  the  inequality 


(3.1) 


l  x  <  1. 


Such  inequalities  are  known  to  define  facets  of  Pj  the  set  packing  polytope 
associated  with  Pj [ 17 ] ;  but  Since  Pj  itself  is  a  face  of  Pj,  it  is  an  open 
question  whether  an  inequality  (3.1)  also  defines  a  facet  of  Pj.  In  this 
section  we  answer  this  question  exhaustively. 

First,  some  definitions  and  basic  concepts.  For  any  polyhedron  P,  let 

dim  P  denote  the  dimension  of  P  (defined  as  the  dimension  of  the  affine  hull 

of  P,  i.e.  of  the  smallest  subspace  containing  P).  An  inequality  nx  <  n0 

is  said  to  define  a  facet  of  P,  if  it  is  satisfied  by  every  x  e  P  and  the 

polyhedron  p”:  =  (x  t  P|*x  =  7tq}  has  dimension  dim  P-1.  If  tix  =  t0  for  all 

x  e  P,  the  inequality  nx  <  is  said  to  define  an  improper  face  of  P.  In 

this  case  of  course  dim  P71  =  dim  P.  To  show  that  *x  <  does  not  define 

an  improper  face,  it  is  sufficient  to  exhibit  a  point  x  e  P  such  that 

tx  <  t  .  Once  this  is  ascertained  to  be  the  case,  dim  p" <  dim  P-1,  since 
o 

(a)  dim  P  is  the  number  of  variables  in  the  system  defining  P,  minus  the  rank 
of  the  equality  system  of  P  (i.e.  of  the  system  of  linear  equations  satisfied 
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by  all  x  e  P);  and  (b)  the  addition  of  the  equation  ix  =  ttq,  not  implied 
by  the  system  defining  P,  increases  the  rank  of  the  equality  system  by  at 
least  1.  Thus  showing  that  rrx  <  wQ  defines  a  facet  of  P  essentially 
amounts  to  showing  that  the  dimension  of  P1,  known  to  be  bounded  by  dim 
P-1,  is  actually  equal  to  this  bound.  The  most  commonly  used  procedure  for 
doing  this  is  to  exhibit  dim  P  affinely  independent  points  x  e  Pn.  Another 
way  of  doing  it  is  to  show  that  the  addition  of  ux  =  ttq  to  the  constraints 
defining  P  increases  the  rank  of  the  equality  system  of  P  by  exactly  one;  in 
other  words,  that  any  equation  satisfied  by  all  x  e  Pn  is  a  linear 
combination  of  the  equations  in  the  system  defining  P1.  In  this  paper  we 
will  take  the  latter  approach,  and  will  use  it  also  to  establish  the  dimension 
of  Pj  itself.  We  will  implement  this  approach  via  a  technique  similar  to  that 
used  by  Maurras  [16],  as  well  as  by  Cornuejols  and  Pulleyblank  [5],  (see  also 
Cornuejols  and  Thizy  [6]). 

We  first  establish  the  dimension  of  Pp 

Let  P  denote  the  feasible  set  of  the  linear  programming  relaxation  of  Pj, 

1  .e. 

_3 

P  =  {x  e  R  | Ax  =  e,  x  >  0}. 

Lemma  3.1.  The  rank  of  the  system  Ax  =  e  is  3n-2. 

Proof.  The  rank  of  Ax  =  e  is  at  most  3n-2,  since  equation  2n  is  the  sum 
of  the  first  n  equations,  minus  the  sum  of  equations  n+l,...,2n-l;  and 
equation  3n  is  the  sum  of  the  first  n  equations  minus  the  sum  of  equations 
2n+l, . . . ,3n-2.  On  the  other  hand,  the  rank  of  Ax  =  e  is  at  least  3n-2,  since 
we  can  exhibit  3n-2  affinely  independent  columns  of  A.  Consider  the  three 
sets  of  columns  indexed  by  the  following  triplets: 
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(2*1,1) »  (3,1,1) . (n.1,1); 

(1.2.1) .  (1,3,1),  ....  (l.n.l); 

(1.1.1) ,  (1.1.2),  ....  (1,1, n) 

The  first  two  sets  contain  n-1  columns  each,  the  last  one  contains  n 
columns.  The  matrix  formed  by  these  columns  (in  the  order  of  their  listing), 
after  deletion  of  the  first  row  of  set  I  and  the  first  row  of  set  J,  becomes  a 
square  lower  triangular  (hence  nonsingular)  matrix  of  order  3n-2,  with  each 
diagonal  element  equal  to  1  and  all  elements  above  the  diagonal  equal  to  0.  j| 
Corollary  3.2.  dim  P  =  n3  -  3n  +  2. 

Proof.  The  dimension  of  P  is  the  number  of  variables  in  its  defining 
system  (n^),  minus  the  rank  of  its  equality  system  Ax  =  e  (3n-2).  || 

We  are  interested  in  dim  Pj.  Since  Pj  c  P,  dim  Pj  <  n^-3n+2,  and 
strict  inequality  holds  if  and  only  if  there  exists  an  equation  ax  =  aQ 
satisfied  by  all  x  e  P;,  that  is  not  implied  by  (not  a  linear  combination 
of)  the  equations  Ax  =  e.  We  will  show  that  no  such  equation  exists. 

n3 

Theorem  3.3.  Suppose  every  x  e  Pj  satisfies  ax  =  ctQ  for  some  a  e  R  , 

a„  c  R.  Then  there  exist  scalars  V  i  e  I,  u V  j  e  J,  and 

o  1  J 

vkt  V  k  e  K,  such  that 

=i  jk  -  vv  v  s  1  * J  *  K 

and 


ao  =  I  xi  +  l  uj  +  I 

0  i  1 1  je J  3  keK 


k" 


Proof.  Let  x  e  Pr«  x  e  {0,1} 
n3 

and  let  x1  e  {0,1}  be  defined  by 


be  such  that  x  .  =  x.  ,  .  =  1, 

1 1J1K1  3J2  2 


xWi 


=  x 1  .  =1 

i  i  k  1  ’ 
1 J2  2 


x*.  .  =  0 

1 lJlkl  13J2k2 


zo 


and  x l =  x . j ^  for  all  other  (i,j,k).  It  is  easy  to  see  that  x  e  Pj 

implies  x'  e  P j ,  and  Since  ax  =  ax'(=aQ)  by  assumption,  we  have 

(3-2)  a.  •  ,  +  ai  -  .  =  a.  .  .  +  a.  .  . 

*  1  rl  1 3^2z  13j1k1  1J2K2 

3 

Next,  let  x  £  PT,  x  e  {0,1}  be  such  that  x.  -  .  =  x.  .  .  =1,  and 

3  1  12J1K1  '  3j2k2 


let  x‘  £  {0,1}  be  defined  by 


*i  Ik  =  *i  ik  * 
3J11  2J2  2 


x-  k  =  xi  ,  k  =  0, 

12J1  1  1  3j2k2 


and  x|jk  =  x1  .k  for  all  other  (i,j,k).  Again,  x  £  Pj  implies  x'  £  Pj, 
and  since  ax  =  ax',  we  have 


(3.3)  a.  •  .  +  a-  .•  k  =  Cl.  .  ,  +  Qt.  .  k  . 

Vri  1  3^2*2  Vri  2J2*2 


Subtracting  (3.3)  from  (3.2)  gives 


(3.4)  a.  •  .  -  a.  •  k 

*  ill  n2JlKl 


=  a  •  •  t  -  a 

1lJ2k2 


i  2^2k2 


If  in  the  definition  of  x,  x1  and  of  x,  x1  we  replace  the  pairs  ( J 2 » k 2 ) 
by  arbitrary  (j,k)  ■  J  x  K,  we  obtain  (3.4)  with  ( J 2 » k 2 )  on  the  right  hand 
side  replaced  by  (j,k).  Thus  the  meaning  of  (3.4)  is  that  two  components  of 
a  whose  index-triplets  differ  only  in  their  first  element,  differ  by  an 


amount  depending  only  on  the  values  of  those  first  elements,  i.e. 


By  symmetry  we  have  equations  analogous  to  (3.5)  for  components  of  a 
whose  index-triplets  differ  only  in  their  second  or  only  in  their  third 


element. 

Now  let  us  denote 

3 i  1  i 2  =  “ijjk  '  ai2J'k 
(3'6)  ’Jjjj  ’  “Ijjk  -  “ij2k 
lk2  "  “ijkj  "  aiJk2 

Then  we  have 


ai jk  '  c,1jj1k1  =  aiJk  "  ai1J‘k  +  “ijjlc  ‘  “ijjjk  +  ai1J1k  ’  01  i i J x k x 

JJL  kkf 
or 


(3.7)  «jjk 

+  v  .  . 

TJJX 

+  5kk  +  ai  i  k  1 

.  V ( i , j , k )  e  I  x  J  x  K. 

Next  we  fix 

Ujl 

and  V. 

K1 

arbitrarily,  we 

let  x.  =  a,  -  k  -  u- 

\  1iJri  Ji 

and  define 

xi 

v 

S"T  ' 

i  £  I  \  { i  j  I , 

(3.8)  UJ  - 

v 

Yjj1  ’ 

j  £  J'ljj}, 

<Jk  = 

v 

5kkL  * 

k  z  K\ { k ^ > . 

Then 

(3.9)  ajJk  .  »,  ♦  Uj  *  vk.  V(I.J.k)  e  I  X  J  X  K . 
Further,  let  x  be  defined  by 

C 1  i  f  i  =  j  =  k 
^  /  0  otherwise. 


( 
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ft 


I 

rt 

< 


« 


* 

9 


l 


l 


< 


Then  k  c  P j ,  *ance  =  a  ,  or 


l  xi  +  L  uj  +  L 
lei  jeJ  J  keK 


V 


Corollary  3.4.  dim  Pj  =  -  3n  +  2. 

Proof.  From  Theorem  3.3,  the  smallest  affine  subspace  containing  Pj  is 
the  one  defined  by  the  system  Ax  =  e;  the  dimension  of  Pj  is  therefore  the 
same  as  that  of  P.  |  | 

Next  we  turn  to  the  constraints  defining  P  and  ask  the  question,  which 
ones  among  these  define  facets  of  Pj. 

Theorem  3.5.  Every  inequality  xs  >  0  for  some  s  £  S  defines  a  facet  of 

PI- 

Proof.  The  statement  is  true  if  and  only  if  the  polytope 
Pj=  (x  e  Pj|x  =  0}  has  dimension  dim  Pj  -  1  =  n3  -  3n  +  1.  Clearly, 
dim  P^  <  n3  -  1  -  r,  where  r  is  the  rank  of  the  system  Asx  =  e,  and  As  is 
the  matrix  obtained  from  A  by  removing  the  column  as.  The  rank  of  As  is 
easily  seen  to  be  the  same  as  the  rank  of  A,  i.e.  r  =  3n-2.  This  is  immediate 
in  the  case  when  as  is  not  among  those  columns  used  in  the  proof  of  Lemma  3.1, 
and  follows  oy  symmetry  for  the  other  case.  Hence  the  dimension  of  Pj  is  at 
most  n^  -  3n  +  1.  To  show  that  this  bound  is  actually  attained,  one  can  use 
the  same  argument  as  in  the  proof  of  Theorem  3.3  to  show  that  any  equation 
ax  =  aQ  (other  than  xs  =  0)  satisfied  by  every  x  £  Pj  is  a  linear 
combination  of  the  equations  Asx  =  e.  The  argument  goes  through  essentially 
unchanged.  || 

As  to  the  equations  of  the  system  Ax  =  e,  each  of  them  defines  an 
improper  face  of  Pj  (i.e.  is  satisfied  by  every  x  £  Pj). 


ft 


ft 


ft 


ft 


ft 


ft 


| 


ft 


ft 
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We  now  turn  to  the  inequalities  (3.1)  defined  by  the  cliques  of  Ga. 


Each  clique  of  class  1  induces  an  inequality  whose  left  hand  sice 
coefficient  vector  is  one  of  the  rows  of  A.  Hence  each  such  inequality  is 
satisfied  with  equality  by  every  x  e  and  therefore  defines  an  improper 
face  of  Pj. 

Next  we  consider  the  inequalities  (3.1)  induced  by  the  cliques  of  class 
2.  Each  clique  in  this  class  is  defined  relative  to  some  index  (triplet) 
s  e  S,  and  has  a  node  set  of  the  form  {s}  U  T(s)  (see  Proposition  2.2).  Let 
C(s)  denote  the  node  set  of  the  clique  of  class  2  defined  relative  to  s. 
Theorem  3.6.  For  n  >  3,  the  inequality 

(3.10)  y  x.  <  1 
teC(s) 

defines  a  facet  of  Pj  for  every  s  e  S. 

Proof .  Let  P,^  =  {x  e  PT|  y  x.  =  1}  and  let  n  >  3.  We  will  show 
C(s)  tcC(s)  1 

that  dim  P^  '  j-  P^.  Let  s  =  00,j0,k0);  then  every  triplet  in  C(s)\{s] 
contains  two  of  the  three  elements  of  s.  Hence  any  x  e  Pj  whose  support 
includes  (iQ,j  ,k  ),  (i  ,j0,k  )  and  (i  ,j  ,kQ),  for  any  i ' ,  i"  t  iQt  j',  j"  f 


jn  and  k  ,  k  f  kn,  satisfies 


)  x 


teC(s) 


^  -  0  <  1.  Thus  (3.9)  does  not  define 


an  improper  face  of  Pj,  i.e.,  dim  P^5^  <  Ppl. 


To  show  that  this  last  inequality  holds  as  equality,  we  prove  that  an 


arbitrary  equation  ax  =  ag  satisfied  by  all  xsPj 
of  the  equations  of  the  system  Ax  =  e  and 


is  a  linear  combination 


(3.11) 


t.C(s) 


xt  =  1. 


Let  x ;  C  ( s ) ,  x  -;  ( 0 , 1 }  c  Rn  ,  be  such  that  x.  .  ,  =  x 

Vri 


1  3  J  2  <  2 


=  1,  where 


( i  1 » j  ^ « k  i ) ,  ( ■>  3  J  2 ' k  2  ^  ^  ~(5)-  Since  n  >  3,  there  exist  at  least  3  disjoint 

triplets  (;,j,k)  such  that  x;^  =  1,  and  at  most  one  of  these  can  belong  to 


i 

1 


24 


C(s);  hence  the  vector  xEP^s)  postulated  here  exists.  Note  that  x  is  of  the 
same  form  as  the  vector  x  used  in  the  proof  of  Theorem  3.3,  and  nence  the  same 
reasoning  can  be  used  to  derive  the  relation 


ai  L  jk 


li,jk 


=  const.,  V(j,k)EJxK:(i1,j,k),(i2,j,k)  i  C(s) 


for  all  i  i ,  i  2  £  I,  analogous  to  (3.5).  Using  the  notation  (3.6),  we  then 
cbtai  n 


ijk  TJJ1  kkL  1 lJlkl 

and  defining  id,  uj,  jeJ,  ,  ksK 

(3.12)  a.jk  =  V •  +  uj  +  vk,  V ( i , j , k )  e 

analogous  to  (3.9). 

Now  let  us  define 


,  V ( i , j , k )  E  S.C(s), 
as  in  (3.8)  we  get  the  relation 
S\C(s) 


ft 


ft 


ft 


l 


I 


.  -  a  •  •  I  —  \  -  U  •  -  Vi 

i  lj  k  i  j  k 
J0  0  J0  0 


(3.13)  =  “i  jk  -  xi  -  "j  -  "k  • 

0  0  0  0 


^k  ai  j  k  "  'i  '  uj  'Jk  ’ 
0  0  0  0 


V  id 


v  jeJ 


V  ksK 


where  (i0,j0,kQ)  =  s. 

Since  teC(s)  if  and  only  if  t  is  either  of  the  form 
(iQ,j,k0),  or  (i0,j0,k),  (3.13)  defines  for  all  (i,j,k) 

w'  11  show  that  i  =  -i^,  V  id,  jeJ,  kiK.  From  (3.13) 

,  .  ,  I  J  K  ,  ,  .  .  I  ,  J 

seen  that  =  -r.  =  ^  ;  hence,  denoting  ,:=  ’ -j  (=  ’  j  = 

0  0  0  '  0  "u 


(i,j0,k0),  or 
£  C ( s ) .  We 
it  is  easily 

-k  >•  n is 
o 


sufficient  to  shew  that 
J  < 

"j  "  Tk  -  ”TC(s)  ’ 


TC  ( 5 ) 


V  id: 


V  jeJ,  k  £  K  . 


it  then  follows  by  symmetry  tnat 


ft 


ft 


ft 


ft 


ft 


ft 
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Let  xeP.^,  xe{0,l}n  ,  be  such  that  x.  .  .  =  x  ,  ,  ,  =  1,  where 

1  oJoo  i  j  k 


I  I  I 


(i  »J  .k  )  i  C(s)  is  chosen  arbitrarily;  and  let  x'e{0,l}  be  defined  by 
x ' ,  =  x*  ,  ,  =1, 

1  J'oko  V  k 


X'.  . 


iik  =  x*  ‘  '  1  "0* 

oJoo  i  j  k 


C(s)  „ 


and  x*  =  x.jk  for  all  other  (i,j,k).  Then  from  xePj  '  it  follows  that 


x'eP^s^,  and  from  ax  =  ax‘(=aQ)  we  obtain 

a.  .  .  +  a  ,  ,  ,=  a  ,  +  at  ,  .  . 

oJo  o  i  j  k  i  jQk0  iQj  k 


Substituting  for  a.  •  ,  a  ,  from  (3.13)  and  a  ,  <,  a  ,  , 

1oJoKo  i  jQk0  i0j  k  i  j  I 


from 


I  I  ^  I  I  I 


(3.12)  (since  i0,j0,k0),(1  ,J0,k0)  e  c(s)  and  O'o'J  »k  )•  O’  ’k  )  *  c(s)) 


we  obtain 

Tt^  +  X  +  U  +  M  +  X1+U|+Vl=TT^,+  \,+  li  +  V  +  X  +  U|+V|  or 

i0  i0  j0  ko  i  j  k  i  i  j0  ko  1Q  j  k 

n-1  =  , ;  and  since  i'  was  chosen  arbitrarily  (subject  to  the  condition 

1o  1 

that  (i  , j '  *k ' )  i  C(s)),  this  proves  that  *  ttJ  =  *C(S).  V  lei. 

J  K  0 

Furthermore,  as  mentioned  above,  *  j  =  TTk  =  itq^,  v  k£K  follows  by 
symmetry. 

We  have  thus  proved  that 


‘ijk 


X.  +  U  j  +  \>k  +  TTc  ( s  )  ’  ^  (1»J»k)  £  C(S) 


xi  +  uj  +  V 


V  ( i ,  j , k )  £  S\C(s) 
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Now  let  X  £  {1,1}  be  defined  by 


1  if  i=iQ  +  e,  j=jQ  +  £,  k=kQ  +  e  (mod  n)  for  e=0,l,...,n-l 
0  otherwise 


Then  X  e  PjS  ,  hence  ax  =  aQ,  i.e., 

“o  =  l  xi  +  l  +  I  v|<  +  ^Cfsl- 

0  iel  1  jeJ  3  keK  K  ^S' 


This  proves  that  the  equation  ax  =  aQ  is  a  linear  combination  of  the 
equations  of  the  system  defining  Pj^;  hence  dim  P^s^  =  Pj  -  1  and  thus 
the  inequality  (3.10)  defines  a  facet  of  Pj.  |[ 

Finally,  we  turn  to  the  inequalities  (3.1)  induced  by  cliques  of  class 
3.  Remember  that  each  clique  in  this  class  is  defined  relative  to  an  ordered 
pair  (s,t)  of  disjoint  triplets,  and  has  a  node  set  of  the  form  {s.t^.tg.tj} , 
where  each  t^,  i  =  1,2,3,  contains  one  element  of  s  and  two  elements  of  t  (see 
Proposition  2.3).  Let  C(s,t)  denote  the  node  set  of  the  clique  of  class  3 
defined  relative  to  the  ordered  pair  (s,t). 

Theorem  3.7  For  n  >  3,  the  inequality 


(3.14)  J  x  <  1 
UeC(s,t)  U 


defines  a  facet  of  Pj  for  all  s,t  e  S. 

Proof.  The  proof  of  Theorem  3.6  can  be  used  to  show  that  dim  < 

Ppl,  and  to  derive  the  relation 


(3.15)  a..^  =  \.  +  uj  +  y ,  V  (i,j,k)  e  S\C(s,t). 


analogous  to  (3.12). 
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Now  let  s  =  (i$,js,ks),  t  =  (it,jt,kt),  and  define 


s  "  Ws  '  \  '  js  '  ks 


(3.16) 


"+  =  “i  i  k  "  xi  "  "  vk 

zl  sJtKt  ’s  Jt  Kt 


~  ai  i  k  ”  Xi  ^k 

2  tJs  t  t  Js  t 


it.  =  a 

z3 


Ws 


'  XU  “ 


To  show  that  it.  =  tt  ,  i  =  1,2,3,  consider  x  e  {0 , 1 } n  , 
li 

that  x.  •  .  =  x.  .  .  =  1,  and  define  x'  by 
sJsKs  'tJtKt 


x£PC^s,t) 


I 


such 


*i  j  k  =  *' 

tJs  s 


i  i  k  ^  * 
sJt  t 


=  xl 


Vsks  1 tJtkt 


=  0, 


and  ^ i jk  =  *ijk  for  a11  other  0»j*k). 


UPCM 


,  and  ax  =  ax'  implies 
+  a-  •  .  =  a.  .  .  +  a. 


oi •  .  *  t  a«  •  i  =  a •  •  i  *r  a .  .  ,  • 

Vsks  Wt  1tJsks  Wt 


Then  x'eP^S,t^  follows  from 


Since  (i^jj^kg),  (ij,j^,k^)  t  C(s,t) 
and  (is,js,kt),(it, js,k$)  i  C(s,t),  we  have  (from  (3.15)  and  (3.16) 

*s  +  xis+  uj$+  vk$+  xit+  ujt+  vkt=  "t^  xi$+  yjt+  vkt+  xit+  uj$+  vk$ 

or  it t  =  it-.  By  symmetry,  n.  =  *=*..,  and  we  are  done .  1 1 
1  t2  s  ^3 

Theorem  3.8  The  inequalities  (3.1)  induced  by  distinct  cliques  define 


distinct  facets  when  n  >  3. 

To  prove  this  theorem,  we  require  the  following  auxiliary  result: 
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Lemma  3.9  For  n  >  3,  there  is  a  feasible  solution  x  with  x$  =  =  xr  =  1 

and  a  pair  of  cliques  with  vertex  sets  C^,  C2,  such  that  1.  s  e 

2.  t  e  3.  r  £  S \ (C2  U  C^);  4.  if  C2  =  C(tQ)  then  t  *  t  and  tQ  n  r  =  0; 

and  if  C£  =  C(t,rQ)  then  r  t  r  . 

Proof.  (i)  Suppose  C]_  =  C(s0),  C2  =  C(tQ)  where 

s0=(is  *ks  to  =  (jt  ,Jt  ,kt  If  to  e  C1  then  w*!«°*9-  let 

”000  000 

1S  =  1t  and  js  =  jt  ,  and  choose  s  =  (isJs  ,k$  )  and  t  *  (i$  tJt,kt  )  with 
00.00  00  ,oo 

i  ^  i  ,  f  j  .  If  t Q  i  then  there  is  an  s  £  such  that 

1  o  0  1 

s  n  t  =  0.  W.l.o.g.  let  s  =  (i  ,j  ,k)  where  k  f  k.  ,  and  choose 

O  3  _  J  _  L 

00  0 

Finally,  select  r  such  that 


s  =  (i_  ,j.  ,kt  )  and  t  =  (if  ,jt  ,k). 
5o  so  0 


t  *Jt 
0  0 


rns=rnt=0.  Then  x  satisfies  conditions  1-4. 

(ii)  Suppose  Cj  =  C(s0),  C2  =  (t0,rQ),  where  tQ  i  Cj_ ,  s0  and  t0  are  as  above 

and  r  =  (1  ,j  ,k  ).  If  rQ  £  Cj  then  let  s  =  rQ,  t  =  t0  and  select  any  r 

0  0  0 

such  that  rns=rnt=0.  Suppose  instead  that  rQ  i  C^.  If 

1  • 

tQ  n  sQ*  0  then  there  is  a  t  e  such  that  t  n  sQ  *  0  and  a 

I  I  I 

corresponding  r  such  that  C2  *  C(t  ,r  );  thus  w.l.o.g.  we  can  assume 

t  n  s  =  0.  Now,  w.l.o.g.,  assume  i  i  i  ,  j  +  j  .  Then  there  is  an 

,  ,  0  0  0  0 

s  e  such  that  s  =  (i  ,j  ,k)  with  k  /  k  ,  k  t  k.  (in  particular, 

5o  so  0  ^o 

if  k  i  k  then  we  can  take  k  =  k  ,  otherwise  there  is  at  least  one 
0  so  5o 

element  of  K  which  is  not  equal  to  kv  or  to  k„  =  k  ,  and  we  can  choose  this 

t  s„  r 

0  0  0 

element  to  be  k).  Choose  s  =  (i$  ,  j$  ,kp  ),  t  =  tQ  and  r  =  (ir  ,jr  ,k) 

000  00 

Then  x  satisfies  conditions  1-4. 

(iii)  Suppose  C1  =  C(s0,wQ),  C2  =  (t0,r0)  where  s0,  t0  and  rQ  are  as  above,  sQ 

i  C2,  tQ  i  C1§  sQn  tQ  »  0  and  wQ  =  (iw  ,  jw  ,kw  ).  Suppose  first  that 

000 

won  lo  =  If  wo  *  ro*  choose  S  =  S0,  t  =  t0.  r  =  W0.  If  W0  =  r0  then 

II  I  I 

C^n  Cg=  0  so  we  can  write  =  C(s  ,w  )  with  s  f  s0,  w  f  rQ.  Then  we 

I  I 

=  s  ,  t  =  t0,  r  =  w  .  Now  suppose  w  n  t  t  0.  It  follows  that 


choose  s 


i  a  -  u-*  j  -  ■  -- w— 


r 
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wQ?<  r  .  Define  a  new  triplet  w  by  selecting  the  element  or  elements  of 
w0\tQ  and  taking  the  remaining  element  or  elements  from  rQ.  By  defintion 
w  n  tQ  =  0. 

I  l  I 

If  w  n  sQ  =  9  and  w  ?  rQ,  choose  s=sQ,  t=tQ,  r=w. 


I  I 

If  w  n  sQ  =  0  and  w  =  rQ  then  choose  s  =  t0,  t  =  sQ, 
r  =  w  ,  exchanging  the  roles  of  and  C2* 


If  was  t  9  then  it  follows  from  the  definition  of  w  that  r  n  s  f  0 

O  0  0 

and  that  the  common  elements  are  in  positions  in  which  wQ  and  t0  share 

common  elements.  W.l.o.g.  suppose  i  ,  *  i  =  i  and  i  =  V  *  Then 

w  ro  so  o  co 

there  is  an  i  e  I  such  that  i  =  i  and  i  }  i  =  i*  .  The  same 

.  o  o  o  o  „ 

holds  for  other  elements  of  w  n  sQ,  if  any.  Define  the  triple  w  from 

w'  by  replacing  i  i  with  i  and  doing  the  same  for  any  other 

w  ro 

•  H 

elements  of  w  n  s_.  Choose  $  =  sn,  t  =  tQ,  r  =  w  . 


For  each  of  the  above  cases,  the  resulting  x  satisfies  conditions  1-4.  || 

Proof  of  Theorem  3.8.  Let  and  Cg  be  distinct  cliques.  Then  t  9 

and  C2\C^  t  0.  From  Lemma  3.9,  there  exists  a  feasible  solution  x  such  that 


y  xs  .  1  and  I  *t  *  1 


SeC 


teC, 


We  will  show  next  how  to  modify  that  solution  to  produce  a  solution  x'  with 


l.  K  *  i 


SeCi 


and  7  x'  =  0 

tic, 1 


I 
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We  construct  x'  by  exchanging  one  element  of  t  and  the  corresponding  element 
of  r,  denoting  the  triplets  thus  created  t',  r'.  It  follows  from  condition  4 
of  the  lemma  that  for  at  least  one  of  the  three  possible  t ' ,  r'  pairs, 

t  I 

t  i  and  r‘  i  ^2*  Since  x^=  Xj  for  all  j  except  x^  =  x^,  =  x^,  =  1 
and  x£  =  x^.  =  0,  x‘  is  feasible  and  satisfies  the  inequality  corresponding 
to  Ci  with  equality  and  the  inequality  corresponding  to  C2  strictly.  Since  x* 
lies  on  the  facet  corresponding  to  Ci  but  not  on  the  facet  corresponding  to 
C2,  it  follows  that  Ci  and  C2  define  distinct  facets  of  Pj.|| 

It  is  easy  to  see  that  each  inequality  x$  >  0  also  defines  a  facet 
distinct  from  those  defined  by  other  trivial  inequalities  or  by  clique- 
inequalities. 


4.  Detecting  Violated  Clique-Facets 

It  is  of  great  interest  in  terms  of  algorithm  development  to  be  able  to 
determine,  for  an  arbitrary  noninteger  solution  to  the  LP-relaxation  of  an 
integer  program,  whether  that  solution  violates  a  facet  of  the  convex  hull  of 
integer  solutions.  One  may  solve  the  LP-relaxation,  then  identify  a  facet¬ 
defining  inequality  that  cuts  off  the  solution  obtained  and  either  add  it  to 
the  constraint  set  of  the  LP,  or  take  it  into  the  objective  function  with  a 
Lagrange  multiplier.  In  general,  for  an  NP-hard  problem  the  facet- 
identification  problem  is  also  NP-hard,  but  for  some  subsets  of  the  facets  it 
may  be  possible  to  efficiently  identify  which,  if  any,  members  of  the  subsets 
are  violated  by  an  LP  solution.  Recent  efforts  to  implement  algorithms  based 
on  this  strategy  (and  employing  branch-and-bound  techniques  when  a  fractional 
solution  is  reached  that  does  not  violate  any  of  the  facets  under 
consideration)  have  met  with  marked  success  [11,  [81.  In  this  section  we 
describe  an  efficient  algorithm  for  detecting  clique-facets  violated  by  an 
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arbitrary  xeP,  i.e.,  an  arbitrary  solution  to  the  LP-relaxation  of  AP3. 
Although  the  cardinality  of  the  set  of  clique-facets  is  0(n°),  (namely,  nJ 
facets  from  cliques  of  class  2,  and  nJ(n-l)J/4  from  cliques  of  class  3),  the 
proposed  algorithm  can  be  shown  to  have  a  worst-case  running  time  of  0(n4). 
In  terms  of  the  number  |S|  of  variables,  this  is  0(|S|4//3). 

We  first  remark  that  given  a  noninteger  x  e  P,  it  can  be  detected  in 

0(n4)  steps  whether  any  inequality  induced  by  a  clique  of  class  2  is 

violated.  Indeed,  each  of  the  n3  cliques  of  class  2  is  associated  with  some  s 
e  S,  and  is  induced  by  a  node  set  of  the  form  {s>  U  T(s),  where  T(s)  is  the 

set  of  those  triplets  that  differ  from  the  the  triplet  s  in  exactly  one 

element.  Since  the  cardinality  of  T(s)  is  3(n-l),  for  each  seS  it  requires 
0(n)  steps  to  identify  and  add  all  x^  such  that  (i,j,k)  e  C(s),  in  order  to 
check  whether  the  sum  exceeds  1  (in  which  case  the  corresponding  inequality  is 
violated)  or  not.  To  execute  this  for  all  seS  therefore  requires  0(n  x  n3)  = 
0(n4)  steps. 

For  cliques  of  class  3  (whose  number  is  0(n®))  the  complexity  bound  is 
not  so  straightforward.  However,  we  will  give  an  algorithm  which  performs 
that  task  too  in  0(n4)  steps.  This  is  possible  due  to  the  following  fact: 
Each  clique  of  class  3  is  of  cardinality  4;  therefore  any  x  e  P  that  violates 
some  inequality  induced  by  a  clique  of  class  3  must  have  at  least  one 
component  of  value  >  1/4.  On  the  other  hand,  we  have 

Lemma  4.1  For  any  x  e  P  and  any  positive  integer  k,  the  number  of 

components  with  value  >  1/k  is  <  kn. 

Proof.  The  value  of  the  linear  program 
(L)  max  (ex  |  x  e  P} 

3 

is  easily  seen  to  be  n,  since  the  vectors  x  e  Rn  and  u  e  R  ,  defined  by 
*s  *  1/n2,  V  s  e  S 
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and 

ur  =  1/3,  V  r  e  R, 

are  feasible  solutions  to  (L)  and  its  dual,  respecti vely ,  with  the  common 
value  of  n;  hence  they  are  optimal. 

Now  if  x  has  more  than  kn  components  with  value  >  1/k,  then  ex  >  n,  a 
contradiction. |  j 

Theorem  4.2.  It  can  be  determined  in  O(n^)  steps  whether  a  given  x  e  P 
violates  a  facet  defining  inequality  induced  by  a  clique  of  class  3. 

Proof.  Let  C(s,t)  be  the  node  set  of  a  clique  of  class  3.  Since 
| C ( s , t ) |  =  4,  if  x  e  P  violates  the  facet-inequality  corresponding  to  C(s,t), 
then  from  Lemma  4.1  x  has  at  least  one  component  >_  1/4.  Further,  if 
C(s,t)  =  {s,t1,t2,t3} ,  from  Proposition  2.4  there  is  no  loss  of  generality  in 
assuming  that  this  happens  for  the  component  indexed  by  s,  i.e.,  that 
x$  >  1/4.  Thus,  instead  of  examining  all  ordered  pairs  (s,t)  such  that  as»at 
=  0,  we  can  restrict  ourselves  to  examine  those  ordered  pairs  (s,t)  such  that 
x$  >  1/4  and  as»at  =  0. 

Consider  now  the  following  algorithm. 

1.  Order  S  according  to  nonincreasing  values  of  x$,  seS. 

2.  For  each  of  the  first  4n  elements  s  =  (is,js,k$)  of  the  ordered  set  S 

such  that  x$  >  1/4  and  each  of  the  (n-1)3  triplets  t  =  (it,jt,kt)  e  S  such 

that  ifc  f  is,  jt  f  j$  and  kt  /  k$,  calculate  the  sum  t(s,t)  =  xi  .  k 

+  x.  .  ,  +  x.  .  ,  +  x.  .  .  .  If  z(s,t)  >  1,  stop:  the  inequality 

sJt  t  Vs  t  tJtKs 

associated  with  (s,t)  is  violated;  otherwise  continue. 

Since  the  algorithm  examines  all  pairs  (s,t)  such  that  as»at  =  0  and 
x$  >  1/4,  it  either  finds  a  pair  whose  corresponding  facet  inequality  is 
violated  by  x,  or  it  stops  with  the  conclusion  that  x  satisfies  all  facet- 
inequalities  induced  by  cliques  of  class  3.  Step  1  is  executed  once  and  it 
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requires  O(n3log  n3)  operations.  Step  2  is  executed  at  most  4n(n-l)3  times, 
and  each  execution  requires  3  additions.  Hence,  the  overall  complexity  of  the 
algorithm  is  0(n4) . | | 
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Intersection  Graph 
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Given  three  disjoint  n-sets  and  the  family,  of  all  weighted  triplets  that  J 
contain  exactly  one  element  of  each  set,  the  3-  index  assignment  (or  3-  dimen- • 
sional  matching)  problem  asks  for  a  minimum-weight  subcollection  of  triplets 
-j  that  covers  exactly  (i.e.,  partitions)  the  union  of  the  three  sets.  Unlike  j 
J  the  common  (2-index)  assignment  problem,  the  3-index  problem  is  NP-complete.  : 
1  In  this  paper  we  examine  the  facial  structure  of  the  3-index  assignment  \ 

.  polvtope  (the  convex  hull  of  feasible  solutions  to  the  problem)  with  the  aid  of 

L_ - - — - ! 


DO 


PC4*** 

t  .an  rj 


1473 


coition  of  1  nov  «*  is  oasouc 
S/ n  0  io  ;-o  i«-  «#o  1 


T* 


I 


I 


» 


| 


» 


» 


SECURITY  CL  ASSlYICATlON  0»  TNIS  RAGE  rfhon  Ooio  e  moron) 


r 


ABSTRACT  (continued) 

the  intersection  graph  of  the  coefficient  matrix  of  the  problem's  constraint  set. 

In  particular,  we  describe  the  cliques  of  the  intersection  graph  as  belonging  to 
three  distinct  classes,  and  show  that  cliques  in  three  of  the  three  classes  induce 
inequalities  that  define  facets  of  our  polvtope.  Furthermore,  we  given  an  ')(n4l 
nrocedure  (note  that  the  number  of  variables  is  n3)  for  finding  a  facet-denning, 
clique- inequality  violated  bv  a  given  noninteger  solution  to  the  linear  programme 
relaxation  of  the  3-index  assignment  problem,  or  showing  that  no  such  inequality 
exists. 
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